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Abstract 

In the Top-down scenarios, the decay of super-heavy particles (to ~ 10 12-16 GeV), situ- 
ated in dark-matter halos not very far from our Galaxy, can explain the ultra-high-energy 
(UHE) cosmic-ray spectrum beyond the Griesen-Zatasepin-Kuzmin cut-off. In case the 
dynamics of this decay is governed by the minimal supersymmetric standard model, a 
major component of the UHE cosmic-ray flux at PeV-EeV energies could be given by 
the lightest neutralino Xii that is the lightest stable supersymmetric particle. Then, the 
signal of UHE x5' s on earth might emerge over the interactions of a comparable neutrino 
component. We compute the event rates for the resonant production of right selectrons 
(eji) and right squarks (qji) in mSUGRA, when UHE neutralinos of energy E x > 10 5 
GeV scatter off electrons and quarks in an earth-based detector like IceCube. When the 
resonant channel dominates in the total Xi e > XiQ scattering cross section, the only model 
parameters affecting the corresponding visible signal rates turn out to be the physical 
masses of the resonant right-scalar and of the lightest neutralino. We compare the ex- 
pected number of supersymmetric events with the rates corresponding to the expected 
Glashow W resonance and to the continuum UHE vN scattering for realistic power-law 
spectra. We find that the event rate in the leptonic selectron channel is particularly 
promising, and can reach a few tens for a one-year exposure in IceCube. Finally, we note 
that UHE neutralinos at much higher energies (up to hundreds ZeV) may produce sneu- 
trino (p) resonances by scattering off relic neutrinos in the Local Group hot dark halo. 
The consequent v burst into hadronic final states could mimic Z-burst events, although 
with quite smaller conversion efficiency. 



1 Introduction 

The presence of a ultra-high-energy (UHE) component in the cosmic-ray spectra beyond the 
Griesen-Zatasepin-Kuzmin cutoff (GZK) [I], as revealed initially by the experiment Fly's Eye 
|2], subsequently by AGASA |Hj, and marginally by Hires |3|, is presently an open problem in 
high-energy astrophysics. The highest-energy particles in the cosmic-ray spectrum are appar- 
ently protons rather than photons, because of their characteristic hadronic-shower fluorescence 
light shapes. Beyond the GZK energies, protons, slowed down by the cosmic Black Body Ra- 
diation, have energy-loss length well below 50 Mpc. This implies that possible sources of 
these extremely-high-energy particles should be within our Galaxy, or inside our Local Group 
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of Galaxies or nearby Clusters. However, there are no such powerful sources in Milky Way 
or nearby Virgo Group, that are correlated with the arrival map of the UHE cosmic rays 
(UHECR). There is a timid, and yet unconfirmed, clustering of AGASA data along the Super- 
Galactic Plane, not sufficient to probe the UHECR confinement in observable GZK volumes. 
On the contrary, UHECR events exhibit isotropy, a signature reminiscent of most distant cos- 
mic edges. It is not clear whether the UHECR spectrum is really showing a GZK cut-off [31 Hj. 
In order to make our UHECR understanding even more confused, there are now first evidences 
of clustered UHECR events by AGASA correlated with known BL Lacs sources 0. This 
UHECR-BL Lacs connection calls for the UHE Blazing Jet of active galactic nuclei (AGN) as 
a source of particle acceleration. Those Bl Lacs sources are mostly at distances above GZK 
lengths. 

Although the GZK puzzle is still experimentally not completely settled, in the past and 
also recently, there have been quite a few suggestions about possible sources and mechanisms 
for accelerating particles to such extreme energies [Hl[7j. 

A class of proposals, generically called Top-down scenarios, tries to explain the high-energy 
part of the spectrum through the decay of super-massive particles (m ~ 10 12_16 GeV) El 
Hm ITU IT2*] , or by their eventual binary collapse and annihilation in place of time-tuned decays 
|13j . These particles are generally predicted in Grand Unified Theories (GUT). Their decay 
lifetime should be comparable with the present age of the Universe, so that their decay products 
would give rise today to the highest energy part of the cosmic-ray spectrum. Also topological 
defect decays can give rise to the UHECR primaries [Bj. In any case, the decay dynamics of 
the super-heavy particles is controlled by the model assumed for particle interactions. In the 
standard model (SM), the decay after fragmentation finally results into photons, neutrinos and 
hadrons. In supersymmetry (SUSY), UHE secondary stable neutralinos (x?) may be born, too 

EH HUES- 

The minimal supersymmetric standard model (MSSM) is a well motivated extension of 
the standard model of elementary particles [HJ. SUSY stabilizes the Higgs-boson mass under 
radiative corrections by introducing an entire spectrum of supersymmetric partners for the 
SM particles. The MSSM predicts a stable weakly-interacting heavy neutral particle, called 
neutralino which can be also a good candidate for cold dark matter 2 . If the masses of 
supersymmetric particles are at the TeV scale (as needed to stabilize the Higgs-boson mass), 
we expect to discover them in a few years at the CERN LHC [T^]. I n the case SUSY is the 
correct theory of particle interactions, it will govern also the decays of super-heavy particles. 
In [TUl rm IT2*] . the super- heavy particle decays in the MSSM extension of the SM have been 
investigated in detail. These decays would involve the entire spectrum of SUSY partners. At 
the end of the decay chain, one would be left with a spectrum (calculable in the MSSM) of all 
the stable particles, that is SM particles plus stable neutralinos. Hence, a crucial prediction of 
these models is the presence of neutralinos in the UHECR spectrum produced by super-heavy 
objects. 

Another possible explanation for the GZK paradox considers relic light neutrinos as a 
calorimeter in resonance with incoming UHE (ZeVs) neutrinos. The latter may come freely 
from distances above GZK from far BL Lacs (hence explaining the UHECR-BL Lacs correla- 
tions) or gamma ray bursts (GRB). They may also arise from topological defects spread all 
over dark halos (Galactic or extra-galactic). This model, the so-called Z-burst model ^B] will 
be discussed in the following in view of a new and analogous sneutrino-bmst model. 

2 In general, in the MSSM, there are 4 neutralinos. Here, we are talking about the lightest one \i- 
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According to the discussion above, we shall assume in this study a rather hard neutralino 
spectrum (i.e. ~ E~ 1,s ) that should be typical for the GUT-relics fragmentation. In order to 
estimate the effects of model dependence in the spectral index, we will also compare our results 
with the ones corresponding to a softer Fermi-like flat spectrum, going as ~ E~ 2 . We shall 
not consider a much harder Z-burst— like spectrum (~ E' 1 ) that is usually used in different 
context, although such a spectrum could be of relevance for a sneutrino-bwcst model. 

One of the major goals for the upcoming UHECR experiments (like Amanda, IceCube, 
Antares, Nestor) is to verify these predictions. In jl7., the possible signatures of up-going 
UHE neutralinos in the future satellite-based detector EUSO [TH] have been investigated. The 
resonant production of quite heavy squarks (m ~ 0.8 — 1.2 TeV) in the xi scattering off 
nucleons in the earth's atmosphere, after crossing all the earth volume, has been considered. 
Heavy squarks are selected so that the neutralino interaction cross sections are moderate, and 
the neutralino flux is not completely stopped by the earth screen (while the corresponding 
neutrino background does not survive the passage through the earth). An analogous process 
leading to a smaller neutrino earth opacity has been also considered earlier, while discussing 
the up-going r neutrinos [TH] . Note that, in EUSO, the UHE neutralino interaction in the 
atmosphere is polluted by more abundant up-going-horizontal r air-showers induced by UHE 
v T earth skimming in the terrestrial crust [2*Uj . 

In the present article, we will focus on the interaction of the highest energy neutralinos 
with both electrons and nucleons in a earth-based detector like IceCube [2T]. In particular, we 
will study the possibility to detect the resonant production of the SUSY partners of electrons 
and quarks, when their masses are larger than the present experimental limits derived from 
direct searches at high-energy colliders (m g > 100 GeV for selectrons [22], m g <^ 300 GeV for 
squarks 3 |2S1)- We will also take into account the constraints on the supersymmetry parameter 
space coming from bounds on the relic density of the lightest neutralino, in case the latter is 
responsible for most of the cold dark matter in the universe. 

Neutralinos with appropriate energy can produce a right/left selectron (e^/ex) resonance 
in the detector via the scattering off the detector electrons 

xi + e~ -> e R)L -> X 

(see Fig. (TJ. Here, X stands for some visible final state arising from the e^L decay, that will 
be eventually revealed by the detector. As far as the scalar mass is not too heavy and xi 
has a relevant gaugino component, for resonant xi energies this channel gives the dominant 
contribution to the scattering cross section, and one can neglect t-channel effects. 
The kinematics and signature is somehow similar to the W~ resonant production in UHE 
antineutrino-electron scattering j21] 

v e + e~ -> W~ -> X , 

that is expected to occur at the Glashow resonant energy E v = M^ v /(2m e ) ~ 6.3 PeV. Since 
we can predict the exact characteristics of this channel (like its hadronic and electromagnetic 
showering and its final visible energies), the corresponding events could be used to calibrate 
the SUSY resonance observation. 

Right /left- squark (g^/g^) resonances can be produced through the neutralino scattering 
off the quarks in the detector nucleons 

Xi + 1 -> Qr,l -> X . 

3 Even lighter squarks are allowed for m§ 3> trig. 
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Although, at the partonic level, this process has the same resonant structure as the selectron 
channel, due to the quark momentum spread inside the nucleon, the cross-section resonance 
peak is smeared, and the Xi incident-energy resonant characteristics are lost. 

The background coming from the UHE v interaction with nucleons [v + q — > X) by t- 
channel processes could create a major continuous background overlapping or overshadowing 
the squark signal. In the following, after computing the relevant Xi event rates, we will address 
the issue of disentangling the SUSY signal from the SM resonant and continuous background 
given by the competitive UHE neutrinos. 

Higher-energy neutralinos [with E x ~ m?/(2m l/ )] could interact with light relic neutrinos 
in the Local Group or Super-Galactic Plane, and produce sneutrino v resonances in the process 

I2S] 

The following sneutrino decay into final states containing hadrons can in principle mimic a 
Z-burst event [HI]- in this paper, we will also go through a v-burst scenario. 

While calculating the neutralino scattering cross-sections, we will stick to a particular 
model of SUSY breaking, namely the minimal supergravity mediated model (mSUGRA) |2fij . 
The mSUGRA and its variants are theoretically well motivated. The parameter space of 
mSUGRA is characterized by five quantities: m and m^, the universal scalar and gaugino 
masses, respectively; A , the universal trilinear scalar coupling (all taken at the scale of grand 
unification); tan/5, the ratio of the vacuum expectation values of the two Higgs doublets; the 
sign of [i, the Higgs mixing parameter (W 3 11H1H2, where Hi, H 2 are the two higgs superfields 
in the superpotential W). While the sign of fi is arbitrary, its magnitude is determined by 
requiring the radiative electroweak symmetry breaking. 

In most of the mSUGRA parameter space, the lightest neutralino is almost a pure bino 
B, the superpartner of the U(l) gauge field 4 . This is particularly important for Xi being 
a good candidate for cold dark matter j2Zj- Then, the right slepton and right squark are 
maximally coupled to the lightest neutralino, and resonant Xif ~~ * fn cross sections are in 
general enhanced with respect to Xif fh cross sections (where / = e,q). 

We will then restrict our study to the right-scalar resonant production and decay into 
the same initial particles Xif- The corresponding production rates will turn out to be quite 
straightforward to compute, the only critical parameters being the physical masses of the 
particle involved in the process (that is the resonance and neutralino masses, nif and m x ). 
In any case, one should keep in mind that the model predicts further new channels associated 
to the left scalar resonances (sort of twin shadows of the right scalars) whose rate is in general 
suppressed with respect to the right scalar by both the lower (and in general more model- 
dependent) relevant decay branching fractions (see Section 3), and the heavier resonance masses 

(since in general m? > m ( - ). 

i L i R n ~ n 

The neutralino energies that are relevant for the resonant reactions XiQ ~^ OS Xi e ~* e i 

x\v — > v are strictly determined by the resonance and target masses [E x ~ mjr/ (2m target )}, 
that is 

9 

Tfl ~ 



4 In general, a neutralino is a linear combination of the two Higgs superpartners Hi, H2 (higgsinos), a W3 
(wino) and a B (bino). 
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E x ~ -^~10 7 -10 9 GeV, (2) 
x 2m e K 1 

9 

E y ~ ~ 10 14 GeV, (3) 

respectively (m p ,m e ,m u being the proton, electron, relic-neutrino masses), where we assumed 
m~ q > 300 GeV, m- e ~ 100 - 1000 GeV, m 9 ~ 200 GeV, and m„ ~ 0.1 eV. 

The plan of this paper is the following. In Section 2, we discuss our assumptions on the 
UHE-neutralino spectrum. In Section 3, we present the cross-sections for the UHE-neutralino 
interaction with electrons and quarks leading to selectron and squark resonances. In Section 4, 
we present the event rates for the selectron and squark resonant production by UHE Xi f° r 
a 1-km 3 ice detector, in benchmark mSUGRA scenarios that takes into account cosmological 
bounds on the neutralino relic density. They will then be compared with the number of 
events expected to arise from the Glashow resonance v e e — *> W~ , and from the v{p) — nucleon 
interaction via charged- and neutral-current processes, assuming comparable UHE-neutrino 
and neutralino spectra. In Section 5, we address the issue of disentangling the SUSY events 
in the detector by looking at their peculiar energy spectrum and ffuence for visible decay 
products. In Section 6, we will discuss a possible v- burst scenario, corresponding to the resonant 
interaction of UHE x°'s with light relic neutrinos in the Local Group or SuperCluster hot dark 
halos. Finally, in Section 7, we present our conclusions. 



2 The UHE-neutralino spectrum 

In this section, we will discuss our assumptions regarding the neutralino flux in the UHECR. 
In our study, we are mainly concerned with models of SUSY-driven decays of super-heavy 
particles. Different studies in the literature show that the final neutralino spectrum is quite 
dependent on the specific fragmentation model jTUHHJUll- Here, we will not adopt any specific 
fragmentation model, but we will follow a more phenomenological approach. We will try to 
keep our study as model independent as possible by assuming for the neutralino flux a —1.5 
spectral index, that should be typical of the fragmentation of heavy objects. The results of 
a switch to a —2 index are also analyzed, in order to estimate the effects of the expected 
uncertainty in the predicted spectral index. Regarding the flux normalization, we will be only 
guided by present and projected experimental limits on UHE neutrino fluxes. 

Neutrinos and neutralinos of comparable energies , while interacting with matter, will pro- 
duce somewhat similar signals in an earth-bound detector. In general, neutralino-interaction 
rates are expected to be smaller than the neutrino ones. This implies less stringent exper- 
imental bounds on neutralino fluxes with respect to neutrino fluxes. Although there would 
be some room for a more abundant incoming UHE-neutralino flux, we will assume, for the 
sake of simplicity and with a conservative attitude, comparable UHE neutralino and neutrino 
fluxes, that is X ~ <f>„. We will then directly turn present neutrino flux bounds into neutralino 
flux limits. Concerning the range of energies we are interested in for selectron and squark 
resonances (i.e., E x ~ (1 — 10 3 ) PeV), AMANDA gives the most stringent experimental limit 
on neutrino fluxes for energies lower than 5 PeV [2H1- in the higher range, upper bounds on 
the UHE neutrino flux comes also from AGASA Q, Flys Eye 0, and RICE [2H|. 

In particular, from the study of a possible excess in cascade events and assuming relative 
flavor abundances u e : : v T — 1 : 1 : 1, AMANDA sets a limit on the total neutrino 
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flux of El g ~ 8.6 x 1(T 7 GeV cm- 2 s" 1 sr" 1 , for 50 TeV< E v < 5 PeV, and ~ 1.5 x 
lO" 6 GeV cm" 2 s" 1 sr" 1 , for 1 PeV< E u < 3 EeV [2H1- 

For comparable interaction strength of neutrinos and neutralinos with matter, neutralino 
fluxes as high as the above limits are allowed, while weaker neutralino interactions allow even 
larger spectra. 

In our study, we will assume that the neutralino flux is comparable with the flux of just 
one neutrino species, i.e. 4> x ~ <j> Vi +vi (* = e il JL i T )- This hypothesis is somehow conservative, 
since, in models of super-heavy particle decays, one can have scenarios where the neutralino 
fluxes are even larger than the sum of all the neutrino components in the upper energy range 

Regarding the neutralino flux shape (or spectral index), we will model it according to the 
typical behavior of present predictions for neutralino fluxes arising from the fragmentation of 
super-heavy objects [TT H ITT | IT2]. 

Models of UHECR fluxes arising from super-heavy X particle decays have been mainly 
motivated by the GZK puzzle. In order to have an appreciable X decay rate today, one has to 
tune the X lifetime to be longer (but not too much longer) than the age of the universe, or else 
store short-lived AT-particles in topological vestiges of early universe phase transitions. Then, 
the estimate of the X lifetime and volume density is very model dependent. Furthermore, the 
internal mechanisms of the decay and the detailed dynamics of the first secondaries do depend 
on the exact nature of the particles. 

Consequently, no firm prediction on the expected flux of neutralinos can be made. However, 
if there are no new mass scales between M SU sy ~ 1 TeV and M x , the squark and sleptons 
behave like their corresponding supersymmetric partners. Then, one can infer the gross features 
of the X-particle cascade from the known evolution of quarks and leptons. 

The upper edge of the neutralino distribution would be then connected to the super-heavy 
particle mass Mx- In the analysis of a super- massive object decay in [T2], the neutralino 
spectrum goes down with the energy as E" 1 ' 3 , and extends up to about E ~ M x /2. 

At large energy fractions x = 2E/M X , the Xi fl ux tends to be comparable or even larger 
that the v flux. In particular, for x > (0.1 — 0.5) (depending on the particular fragmentation 
treatment [TQ1 El E|) neutralino and neutrino fluxes are in general comparable, while at 
x < 0.1 neutralinos are depleted with respect to neutrinos. Also, the total energy carried away 
by the neutralinos can be comparable with that of the neutrinos. In fact, when the initial 
decay of the X particle is into SUSY particles, the total energy of the decay carried by the 
neutralinos can be even larger than the neutrino energy (see, e.g., Tables 2.1, 2.2 and 2.3 in 
the last reference in jJJ]). 

Consequently, the neutralino flux in the energy range relevant for resonant production of 
selectron and squarks [i.e. E x ~ 10 6 ~ 9 GeV, see eq. (HJ),©], could be quite depleted when 
supposing Mx ~ 10 12-16 GeV in order to explain the GZK puzzle, implying neutralino energy 
fractions x < 0.001. 

For instance, assuming Mx — 2 x 10 12 GeV to fit the GZK cosmic rays data (see Fig. 1 
in J7j), and using the fragmentation model in [T2| . one gets neutralino fluxes of the order 
of E\ = 10~ 8 GeV cm" 2 s" 1 sr" 1 at E x ~ 10 2 PeV (see Fig. 2 in [T7j). This flux would 
be more than one order of magnitude smaller than the value that can be extrapolated from 
present experimental upper bounds on neutrino fluxes at lower energies, assuming, as in jTTj, 
a E~ 13 spectral behavior. 

On the other hand, a lighter X-particle mass in the range 10 8-9 GeV (as considered for 
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instance in |3T)j) would provide a more appropriate setting for resonant selectron and squark 
production, enhancing the neutralino flux up to a level close to neutrino fluxes in the 10 7-9 
GeV energy range. A Mx value in the 10 8-9 GeV ballpark would not be in conflict with any 
theoretical or experimental consideration, although it would be less interesting in view of a 
possible explanation of the GZK puzzle. 

Anyhow, we stress that at the present stage one could also think of models the can reconcile 
the 'low' energy range needed for resonant selectron and squark scattering with the 'large' Mx 
value fitting the GZK energy range. For instance, assuming for the heavy X object a lifetime 
shorter than the one of the Universe, one would presently expect a tail of redshifted secondaries. 
If the X decay occurs smoothly at redshift 10 3 > z > 1, then the peak X decays at largest 
redshift (z ~ 10 3 ) even in the case Mx ~ 10 12 GeV would eject its neutralino component at a 
lower energy of the order E x ~ Mj/[2(1 + z)\ ~ 10 9 GeV today. 

All the previous considerations converge in our choice of a "reference" neutralino flux given 
by a fragmentation-like power law with spectral index —1.5, and normalization well compatible 
with the AMANDA bounds : 



El 4^ = 1(T 7 GeV cm' 2 s" 1 sr~ 
x dE v 



E 



x 



10 7 GeV 



0.5 



(4) 



Although one could easily build up models matching the above "reference" normalization, 
in the following analysis it will be straightforward to rescale our results on expected event 
numbers according to a possible more conservative (smaller) normalization factor in eq. (j3J). 

In order to assess the effect of the uncertainty in the spectral index arising from differences 
in modeling the super-heavy object fragmentation, we also considered, for comparison, a Fermi- 
like neutralino flux with spectral index —2 : 

E 2 X ^ = 1(T 7 GeV cm" 2 s" 1 sr" 1 , (5) 

Ci J—J y/ 

whose normalization is also compatible with experimental neutrino limits. We than compared 
event numbers relative to eqs. (J1J and 

Summarizing this section, the details of UHE neutrinos and neutralinos fluxes are quite 
model dependent, and we are in a very preliminary understanding of their possible spectra. 
Hence, we decided to use two phenomenological assumptions for the UHE neutralino spectrum 
differing by the spectral index, and reflecting possible differences in modeling the fragmentation 
of a super-heavy object. In both cases, the flux normalization has been fixed on the basis of 
present UHE neutrino experimental bounds. This hypothesis could be quite conservative, 
since the lower neutralino-nucleon cross section (<j x n Jo^^n, see Section 4) would make any 
possible astrophysical source of UHE neutrinos and neutralinos more transparent to (and hence 
more efficient in producing) the Xi component. 

Assuming a different normalization in eqs. (j3J) and (0) will simply imply the rescaling of 
the event numbers obtained in our final analysis. 

3 Resonant cross sections for UHE neutralino interac- 
tions with matter 

UHE neutralinos with proper energy will interact with the electrons and quarks inside the 
detector, and produce selectron or squark resonances via Xi + ( e_ ) l) ~> (e, q) — ► X (see the 
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Figure 1: Feynman diagram of the process Xi + / ~ / ~ * X, where / can be an electron, a 
quark or a neutrino. 



Feynman diagram in Fig. Q where / is either a lepton or a quark). The symbol X stands for 
a given visible final state arising from the / decay. 

Near the resonance, the cross-section a of the process can be written in terms of the particle 
masses and the fx± center-of-mass (cm.) energy s in the Breit-Wigner approximation 

In the equation above, mj is the slepton/squark mass, is its total decay width, and B(f — > 
/ x5); -^(/ — * ^0 are their decay branching fractions [-B(/ —►...) = r(/ — > . . into the 

final states / Xi anc i X, respectively. The modulus of the cm. 3-momentum of the initial 
particles k can be expressed as 

2 2 

m% — mi 

|k| = X - (7) 

where m x is the neutralino mass. Then, the cm. energy s is related to the incident neutralino 
energy E x in the laboratory by 

s = m 2 x + 2m f E x , (8) 
where mj is the target mass. At the resonance (s ~ m 2 j), the neutralino energy is then 

2 2 

m % — mi 

E peak _J_ X , g x 

x 2m f ' v ; 

which sets the neutralino energy scales relevant to our study: PeVs-EeVs for squarks and 
selectrons, and hundreds of ZeVs for sneutrino bursts. 



The peak cross section a peak = a(s = m" 2 -) can be obtained from eq. (JHJ) in a straightforward 



way 



ni 



2 



2 



* peak = 4 ( -J^-J } B(f - / x?) B(f - X). (10) 



m% \ m } — m% 

I \ j * , 



In the final event-rate estimate, a crucial quantity will be given by the product a pea k Tj, the 
resonance peak cross section times its total decay width (the latter setting the actual energy 
width of the resonance curve), that in general will depend on all the 5 mSUGRA parameters. 
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In this analysis, we will concentrate on some important scenarios where a relatively small 
parameter dependence is left in the cross sections, apart from the leading mj and m x mass 
dependence. In particular, in these (quite general) scenarios, there will be little parameter 
dependence in the / decay branching fractions into both the initial fxi arid final X states, 
these branching fractions being the crucial parameters entering a pea k- These conditions most 
naturally realize in processes where the resonance is maximally coupled to the initial fxi state 
[that is B(f — > f xi) — 1]> an d, at the same time, one considers as decay products just the 
same initial fxi particles, in the process 

XX + /-/-X1 + /- (11) 
Then, the corresponding peak cross section becomes 

^^(_^L_V (12) 

In the same scheme, also the parameter dependence of the total width Tj will be mainly 
restricted to phase-space mass effects. 

A natural framework were to implement the above picture is given by the right selectron 
and right squark resonances in mSUGRA. In mSUGRA, the right selectron cr is in general one 
of the lightest superpartner, and the lightest neutralino is mainly a bino B. For conserved R 
parity, the dominant e# decay channel is Cr — > exl, with a corresponding branching fraction of 
almost 100 %. On the other hand, the left selectron e L is usually more coupled to the second 
lightest neutralino x%. an d to the lightest chargino Xi ■ Hence, when allowed by phase-space, 
the decays ei — > x\ e an d — > xt e tend to be dominant. Therefore, B{ei — > e xi) is i n 
general not very large, and depletes the peak e xi ~ * &l — » X cross section. 
On the other hand, for the right selectron, according to eq. (JTUjl . the peak cross-section is 
maximized when looking to at a x1 e final state in the channel 

X°i + e^e R ^Xi + e, (13) 

being rrte R and m x the only relevant parameters that govern the cross section, as from eq. (|12)1. 

Then, the dominant signature of the resonant interaction of a neutralino with an electron in 
the detector will be an energetic electron carrying away a fair fraction of the initial neutralino 
resonant energy. The cr signature will be totally characterized by an electromagnetic showering, 
while the Glashow W resonance will be mostly proceeding through hadronic channels. 

For the squark sector, that is relevant in the neutralino interactions with the nuclei of 
the detector, a similar discussion holds regarding the branching fractions of the right and 
left squarks into qxl [in particular, one has B{q~R — > q Xi) ~ 1, B(qi — > q xf) ~ 2/3, and 
B(c[l —> q X°) ~ 1/3] j unless the gluino is lighter than the squark. In the latter case, the strong 
interacting decay ~^ 19 is by far dominant, and the effective squark coupling to the initial 
qXi is suppressed [due to the small branching fraction B(q — > q xi) ~ T(g — > q Xi)/T(q Qd)]- 
However, (as noted above) a pea k is not the only relevant quantity for event rates, since it enters 
the latter through the product cr pea kTg R that also involves the total squark width. With the 
onset of the gluino decay channel, the decrease of B(q~R — > q xi) in eq. (fTnj) is compensated 
by the corresponding increase in the total width Tq R . Hence, when including the gluino decay 
q~ R — > qg into the visible final states of the process [hence keeping B(q R — ► X) ~ 1] the final 
event rate will be the same as in the case of the lighter-than-gluino squark. 
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In the following we will concentrate on the resonant production of up and down right 
squarks, ur and d R , in the processes 

Xi + u, d -> u R , d R -> x°i + u, d , (14) 

assuming B{u R ,dR — > X? + w ) ^) — 1- The dominant signature in the detector will then be 
given by a strongly interacting particle initiated shower, carrying some relevant fraction of the 
initial Xi energy. 

4 Expected event rate at IceCube in mSUGRA bench- 
mark scenarios 

In mSUGRA with conserved R parity, the lightest neutralino is stable and provides an ideal 
candidate for explaining the origin of the cold dark matter (CDM). Taking into account the 
constraints on the CDM density obtained by combining the precise measurement of the cosmic 
microwave background by the WMAP experiment [31] and other cosmological data (that give 
0.094 < VtcDuh 2 < 0.129 at 2-a level), one should then require that the relic density of 
the lightest neutralino from primordial universe be compatible with QcDMh 2 . In particular, 
requiring that most of the CDM is made up of lightest neutralinos restricts quite a lot the 
allowed region of the mSUGRA (m ,m 1 / 2 ) parameter plane (see, e.g., j32], and references 
therein) . 

In jsjjj], a set of 13 mSUGRA benchmark scenarios (characterized by a given set of values 
for the mSUGRA parameters m , m.i/ 2 , A , at the GUT scale, and by the values of tan/3 and 
sign /i) have been chosen, respecting the following experimental constraints: non-observation 
of SUSY partners at colliders (especially at LEP) [221 EH Ell j bounds on deviation from the SM 
b — > 57 decay rate; agreement of the lightest-neutralino relic density with 0.094 < QcDAih 2 < 
0.129. 

Of the 13 benchmark points in [33], 5 points (B', C, G', I', L') are in the 'bulk' region (i.e. 
at moderate m , rax/2 values), 4 points (A', D', H', J') are along the co-annihilation 'tail' (at 
moderate m and larger m^, where the neutralino-stau co- annihilation channel considerably 
affects the neutralino relic density), 2 points (K', M') are along rapid-annihilation 'funnel' 
(where both mo and mi/2 can grow large with large tan/?), and 2 points (E', F') are in the 
'focus-point' region, at very large mo- 

In the present analysis, the two crucial parameters ruling the signal rates are the resonance 
scalar mass mj and the splitting in the square of the scalar and neutralino masses, m 2 j — m 2 ,. 
Then, we found that scenarios like K', M' ('funnel') and E', F' ('focus points'), where the 
scalar mass parameter is quite large (mo > 1 TeV =^> mj > 1 TeV), give rise to drastically 
reduced signals (less than 0.01 events/year in Icecube, see below). On the other hand, for 
such large scalar masses, t-channel diagrams in the Xif scattering can become dominant over 
the resonant s-channel, in case neutralino has a non-negligible higgsino component (cf. 'focus- 
point' scenarios). Then, the approximation of restricting the cross section evaluation to the 
resonant channel (on which the present analysis in based) could fail. Hence, in our study, 
we will consider only the 'bulk' and co-annihilation benchmarks, leaving the analysis of the 
'funnel' and 'focus-point' scenarios to a more complete treatment. 

The 9 benchmark points we include in our analysis are reported in Tabled They correspond 
either to the 'bulk' region (B, C, G, I, L) or to the co-annihilation region (A, D, H, J). In all 
scenarios, A = 0. 
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Model 


A 


B 


C 


D 


G 


H 


I 


J 


L 


mi/2 (GeV) 


600 


280 


400 


525 


375 


935 


350 


750 


450 


m (GeV) 


107 


57 


77 


101 


110 


233 


180 


298 


303 


tan f3 


5 


10 


10 


10 


20 


20 


35 


35 


47 


sign(/i) 


+ 


+ 


+ 


— 


+ 


+ 


+ 


+ 


+ 


tth z 


0.128 


0.123 


0.122 


0.116 


0.110 


0.123 


0.117 


0.119 


0.113 


m x (GeV) 


243 


107 


158 


212 


148 


388 


138 


309 


181 


rrie R (GeV) 


254 


128 


175 


226 


184 


423 


227 


412 


349 


T Sr (MeV) 


9.13 


57.2 


29.7 


16.3 


103 


53.4 


A A O 

448 


396 


931 




21.1 


6.58 


9.34 


13.3 


2.57 


2.17 


0.477 


0.301 


0.150 


Men 


28 


32 


26 


23 


8.5 


2.3 


1.9 


0.55 


0.42 


(km yr ) 


37 


45 


34 


30 


7.7 


1.4 


1.1 


0.20 


0.14 


mq R (GeV) 


1194 


596 


825 


1057 


781 


1798 


748 


1487 


961 


r 9 - fl (GeV) 


2.45 


1.23 


1.70 


2.17 


1.61 


3.65 


1.55 


3.04 


1.99 


< afc (nb) 


7.43 


28.6 


15.3 


9.46 


17.1 


3.32 


18.5 


4.81 


11.3 




0.35 


3.0 


1.1 


0.52 


1.3 


0.08 


1.5 


0.17 


0.70 


(km 3 yr x ) 


0.22 


2.6 


0.89 


0.36 


1.1 


0.04 


1.2 


0.09 


0.51 



Table 1: Definition of different mSUGRA scenarios and corresponding event-number expecta- 
tions Me R ,q R i for resonant neutralino scattering in ice. In all scenarios, we assume A = 0. The 
resonance decay widths and peak cross sections for the right selectron and right squark are 
also shown. The upper of the two entries in the Me R ,q R raws corresponds to the event number 
calculated by the neutralino flux in eq. (JU) (/5 = 1.5), while the lower entry refers to eq. (jSJ) 
(flo = 2). When computing the squark event rates, an incident neutralino energy threshold of 
1 PeV is assumed. The relic DM density flh 2 corresponding to each scenario is also shown. 

Starting from the values of m , m 1 / 2 , A at the GUT scale, and tan/3, signfi, all the 
physical masses and couplings of the spectrum of SUSY partners can be calculated at the 
electroweak scale by solving a set of renormalization group equations. There are many public 
spectrum calculators doing this task, including second-order radiative corrections (see, i.e., 
[3*5]). The differences in the corresponding results is considered a good estimate of the present 
theoretical uncertainty in the computation of SUSY spectra. 

Then, the high precision of the WMAP QcDAdh 2 determination (that slimmed considerably the 
allowed mSUGRA parameter regions) gives rise, through the prediction of the neutralino relic 
density versus mSUGRA parameters, to some sensitivity of the allowed mSUGRA regions to 
the choice of the spectrum calculator in the analysis. The differences in results is considered a 
good estimate of the present theoretical uncertainty of this kind of analysis. 
The impact of the sensitivity to the details of the mSUGRA spectrum has been analyzed in 
[3*o] through the code micrOMEGAs 1.3 [HHl that can evaluate the neutralino relic density 
assuming different codes for the mSUGRA spectrum. One finds that mass differences of about 
1% in the mSUGRA spectrum can imply sometimes a 10% variation (and even larger for high 
tan /3 and mo) in the corresponding neutralino relic density. 

The above discussion motivates the choice of our benchmark points in Table ^ In our 
analysis, we used Isajet 7.71 [3JJ to compute the mSUGRA spectrum 5 , and micrOMEGAs 1.3 
to compute the corresponding neutralino relic density (as implemented in the last reference in 

5 In Isajet7.71, we set m h (m b ) MS = 4.214 GeV, a s (M z ) MS = 0.1172 and the top pole mass at m t = 175 
GeV. 
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Figure 2: Cross sections (in nb) for the resonant right-selectron/right-squark production by 
UHE neutralino scattering on electrons/nucleons, in the four mSUGRA scenarios B, G, I, L 
defined in Tabled For comparison, the dashed lines show the scattering cross sections of UHE 
anti- neutrinos on electrons (when marked by W) and nucleons (when marked by uN). 



[33]). The small differences in the definition of our benchmarks A, B, C, ... in Tabled with 
respect to the corresponding post-WMAP benchmarks A', B', C, ... defined in Table 2 of [33] 
(which are based on ISASUGRA 7.67) are due to the small adjustments in the choice of mi/ 2 
and m needed to pass both the relic-density and LEP constraints by assuming a different 
code. 

In Table one can read the neutralino relic density, Qh 2 , corresponding to each scenario, 
that is well compatible with the assumption that the lightest neutralino be the main CDM 
component. For each benchmark, we show the lightest-neutralino, right-selectron and right- 
squark masses (masses for squark of different flavors u,d,s,c are assumed to be degenerate 
and equal to m^). In all scenarios, one finds < 0.5 TeV, and rriq R < m g , so that 
B(qr ~^ 1 Xi) ~ 1- Total right-selectron and right-squark widths, and Tg R are also shown 
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in Table [TJ along with the respective peak cross sections, a^ e R k and cr^ k , as defined in eq. (fTUjl 
with B(f - X) ~ B(f -/*?). 

In Fig. |2J we show (for the B, G, I, L scenarios) the cross sections for the resonant production 
of a right selectron and a right squark from the UHE neutralinos scattering on electrons and 
nucleons, respectively, versus the incident neutralino energy. 

For comparison, we also present the resonant neutrino cross section for v e e~ — > W~ — *> all, 
and the (non- resonant) neutrino-nucleon cross section. Indeed, a Breit-Wigner W peak cross 
section is expected at the incident neutrino energy (E v ) 

EP eak = M^-ml ^ Myv_ ^ Q 3x 1Q 6 GeV (15) 

although present observation of UHE energy neutrino events does not reach such energies yet. 

Cross sections for the resonant squark production do not show the Breit-Wigner structure. 
Because of the continuous quark momentum distribution inside a nucleon N, any value of the 

incident neutralino energy larger than a threshold [E x > ^ m ^ x , where is the nucleon 
mass) can produce a squark at resonance. In this case, the cross section evaluation involves a 
convolution of the partonic cross section in eq. © (with / = u,d and sea quarks) with a parton 
distribution function f q (x, Q 2 ) , where x is the momentum fraction of the proton carried by the 
parton, and Q is the energy scale at which we evaluate this distribution. In our computation, 
we used the CTEQ4L set of structure functions in [3~8j . 

The ur and da decay widths are small compared to their masses, and one can substitute 
the Breit-Wigner propagator in the partonic cross section by a 5-function. Then, the squark- 
resonance cross section in the x^-nudeon scattering becomes: 



peak jn 

qR {E x ) = tt £ x q f g (x q ,ml) (16) 

where, x q = m $ a /( m x o + 2 i? x mjv). The summation in eq. (jlfij) is over the u- and d-type 
quarks (thus including also s and c quarks) in an isoscalar target. In our scenarios, one has 
m u R - m d R - m c R - m 's R and ~ T &R ~ 4 T iR ~ 4 F Sr . 

The actual event rates Afz can be calculated from the cross section by convoluting it with 
the appropriate neutralino flux, and multiplying the result by the number of the target elec- 
trons/nucleons A/r in a given detector volume, and by the exposure time At, as follows 



r dN 

A/} = M T At j aj{E x ) — dE x . (17) 

X 



For selectrons and squarks aj(E x ) will be given by eq. (0) and eq. respectively. We 

assumed the two different neutralino fluxes in eq. and eq. (J5J). Nt is given by the product 
of the detector mass (in g units) and Avogadro number (Na) for nucleons (times 11/18 for 
electrons) . 

The selectron width being small compared to its mass (see Table [Q, we approximated the 
Breit-Wigner propagator in eq. © by a 5-function. The corresponding expected event number 
is then 

A/L = o Sl , a?T k IV, ( ^ 1 (Er k Y"° , (18) 



'e R ~e R » eR , eR j y^ x j 

where the factor ae R (depending on the mass of the detector, exposure time and nature of the 
neutralino flux) includes the product of At, A/r, and the normalization of the fluxes defined in 
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eq. (jlj) and eq. (JSJ). The exponent (3q sets the shape of the neutralino flux, and is equal to 1.5 
and 2, respectively, for the two cases, while E^ eak , expressed in GeV, is given by eq. © where 

For an ice detector of 1 km 3 and one year of exposure, one has ag R = 7.24 x 10 3 [2.29 x 10 7 ] 
nb _1 GeV^ 0-1 for the flux in eq. (J3J) [eq. (j5|) ]. The large difference in the a factors is an artifact 
due to the exponent /3 in the formula above. 



Number of events from UHE u, v interaction 


W~ -> ev e 


W — > hadrons 


+ Vi)cC+NC 


2 


12 


55 


2.5 


15 


56 



Table 2: Number of events from UHE neutrino/ antineutrino CC and NC interaction in IceCube 
per year. In each entry, the upper number refers to the flux in eq. (/3q = 1.5) for the Ui + Pj 
spectrum, while the lower number assumes the flux in eq. (0) (Po = 2). An incident-neutrino 
energy threshold of 1 PeV is used. 

For the g^-resonances, due to the energy distribution of the partons in a nucleon, the 
number of events will get contributions from a continuous range of incident neutralino energies 

Kn = ^ [ . a, R (E x ) dE x , (19) 

where, for a 1-km 3 ice detector and one year of exposure, aq R = 3.77 x 10 3 [1.19 x 10 7 ] nb _1 
Q e yA)-i f or the flux in eq. (jlj) [eq. (jSJ], and aq R (E x ) is given by eq. (JEJ). 

One can note that the peak cross section in eq. (|T2*|) is enhanced for small mass differences 
between the selectron/squark and the neutralino. However, in the product <j^ eafe that enters 

the event numbers in eq. (J18|) . and in eq. (J19|) through eq. (|16|) . the l/(wr~ — m 2 x ) 2 behavior in 

(j p - eak is compensated by the resonance-width, that goes as (m| — rn x ) 2 - On the other hand, a 
residual enhancement associated to the resonance-neutralino mass degeneracy arises from the 
factor (-E£ 6afc ) ° in eq. (JTBJ) (cf. eq. Q), and from the lower limit E x in = (m| - m 2 x )/(2m N ) 
of the integral in eq. (|19j). Hence, selectrons/squarks and neutralinos that are closer in mass 
tend to increase the expected event number. 

In Tabled we present the number of events {Me R , q R ), per km 3 per year, expected in an ice 
detector for the different mSUGRA benchmarks. The number Mq R was always worked out by 
assuming as lower limit in the integration in eq. (fT9"|) the value E x m ~ 10 6 GeV in order to 
cut off the atmospheric neutrino pollution. In Table ^ the raws corresponding to the number 
of events have two entries. The upper one corresponds to the neutralino flux defined in eq. (Jlj) 
(/5o = 1-5), while the lower one assumes the flux in eq. (JHJ) (/3o = 2). 

For comparison, the expected event rates for the resonant W production from v e e scattering, 
in the hadronic and leptonic channels, are presented in Table 121 Since the fluxes in eqs. (0) 
and (J3J are relative to one flavor of v + u, assuming a complete neutrino flavor mixing due to 
non-zero neutrino masses, in Tabled we use a v e flux that is half of those defined in eq.(J3J) and 
eq.(J5J). We also show total v-N and i>-N charged-current (CC) plus neutral-current (NC) event 
rates. 

In scenarios A, B, C, D, where one has both a light selectron (m e - B < 250 GeV) and quite 
degenerate me R /m x values, one obtains in IceCube a few tens of events per year in the electronic 
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channel. Note that the selectron/neutralino degeneracy favors Me R in the harder-spectrum 
case with /3 = 1-5. In scenarios H, I, J, L the larger splitting in m e - fl and m x , and the quite 
heavy selectron deplete the event rates by more then one order of magnitude (0.2-2 events 
per year). Note that in this case the /3q = 2 spectrum gives higher Me R - The benchmark G 
shows an intermediate case, where the selectron is quite light, while mg H and m x are not much 
degenerate, giving rise to about 8 events per year. 

Concerning the hadronic channel, due to the large squark masses (m^ > 600 GeV) and 
their large splitting with m x , the signal is modest in all benchmarks, reaching the level of 1-3 
events per year only in the B, C, G, I points. 

With the assumed neutralino fluxes, one could expect a clear SUSY signal in IceCube in 
the electronic channel for most of the analyzed benchmarks. Only scenarios J and L fore- 
see less than 1 event per year in both the electronic and hadronic channel. We recall that 
the normalization of the fluxes in eqs. (@J) and (0) is quite conservative. Due to the lower 
neutralino- nucleon cross section a xN < j^a v ^ (cf. Fig. |2J), one is presently allowed to assume 
flux normalizations that are an order of magnitude larger than the ones adopted here with- 
out violating any experimental constraint. This would make the event statistics even more 
promising. 

5 Disentangling the neutralino signal by neutrino cali- 
bration 

We will now briefly address the problem of the effective visible- energy spectrum coming from 
the two body decays fn — > fxi, that is of course a crucial characteristic in the detection of a 
resonance signal. 

In the resonant e# decay into an electron plus an invisible neutralino Xi> only the electron 
energy contributes to the visible energy. In Fig. El we plot the visible-energy spectrum Az^j 
(where z = \og w E vis , and Az = 0.1 assumes a 10% experimental resolution) corresponding to 
the resonant e# production and decay into ex?, for benchmarks A, B, G, I. The positions of the 
sharp end-points of the spectra are determined by the resonance and neutralino masses. For 
comparison, we also show the expected visible-energy spectrum for a resonant W decaying into 
either ev or hadrons. When the W decays hadronically, the total energy of the initial incident 
neutrino gives rise to visible energy in the final state. In this case, the Az^j spectrum has a 
resonant structure corresponding to the initial neutrino energy E v ~ m^/(2m c ) ~ 6.3 x 10 6 
GeV. When the W decays leptonically, only a fraction of the incident neutrino energy goes into 
visible energy, whose distribution is characteristic of a two-body decay, just as in e# — > e\\- 
The different positions of the spectrum end-points could then provide a handle to disentangle 
the visible energy of the e# decay from the one of the W decay. 

The hadronic W~ decay defines a sharply clustered group of events [about 12 for the 
Po = 1.5 spectrum, and 15 for the (3q = 2 power law, per year, cf. Table EJ, whose intensity 
may be correlated with the two electromagnetic channels arising from W~ — > e~ + u e [2 and 2.5 
events] and W~ — > t~ + v T [2 and 2.5 events]. The channel W~ — > + z/ M [2 and 2.5 events] 
is delivering just a single PeV muon track, that is quite difficult to recognize. The electron 
and tau electromagnetic energy spectra grow linearly with energy up to E = 6.3 PeV, and 
their average energy is about half of the hadronic-channel energy. The distribution is spread 
according to a continuous power law, as shown in Fig. El The r birth and decay will be in 
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general source of a characteristic double-bang event |3U] . 

In conclusion, the appearance of a sharp edge at some energy different from the expected 
W~ hadronic, electronic, and tau bumps could point to some new physics corresponding to the 
resonant e~R production from UHECR neutralinos. Such reconstruction seems possible for a 
neutralino flux that is comparable to (and, in favorable scenarios, even lower than) the present 
lower bounds on neutrino fluxes. 

The mass ranges considered are of great astrophysical interest, and within the reach of the 
CERN LHC. A near-future underground i/-detector could then cross-check a possible SUSY 
discovery at future colliders. 



l 1 1 1 

W— >hadrons 




Log 10 [E vis (GeV)] 

Figure 3: Visible energy spectra of the resonance decay products [assuming the flux in eq. (jlj) 
(/?o = 1-5)]. For the resonant W production and decay, the dash-dotted (dashed) histogram 
represents the electron (hadronic) energy spectrum in the lab. The four solid-line distributions 
refer to the electron energy spectra arising from e# — > ex° in the mSUGRA scenarios A, B, G, 
I defined in Table [3 
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6 The role of the sneutrino burst 



The resonant production of Z-bosons via UHE v scattering off massive relic neutrinos has 
been the backbone of the so called Z-burst mechanism ^B] that could explain the GZK 
puzzle. UHE neutrinos in Z-burst are the messengers able to overcome the GZK cut-off, 
by offering a link between UHECR and most distant cosmic sources. About 70% of the 
Z— bosons decays hadronically, thus producing protons, neutrons and anti-nucleons as well 
as pions (and consequent photons from neutral-pions decays). The Z-decay products could 
have quite naturally energies beyond the GZK scale for appropriate relic-neutrino masses, 
since E v ~ Mf/(2m„) ~ 4 • 10 22 ( °^ eV ) eV, which corresponds to a UHECR proton energy 
of E p ~ ^ • Mz/(2m u ). If relic-neutrino halos are present at distances less than about 50 
Mpc 6 , then these protons or anti-protons (or photons) could arrive to the earth atmosphere 
without much energy loss. Consequently, air-showers produced by these extremely high-energy 
nucleons could give rise to the super GZK events, as detected by the AGASA experiment |3] • 
Crucial ingredients of this picture are of course sufficiently large UHE-neutrino fluxes and relic- 
neutrino densities in the Local Group or Super Galactic Plane. Then, the Z-burst model has 
a good potential for explaining observed correlation of UHECR arrivals from BL-Lac sources 
at distances well above GZK, and electron pair losses [5]. 

In this section, we will consider the resonant sneutrino (u) production from the UHE 
neutralinos scattering off relic neutrinos, and its subsequent decay into visible final states 

X°i" -> v -> X vis , (20) 

where X vis stands for the all final states with at least one visible particle. In particular, we will 
discuss whether this process can have a role similar to the Z-burst in the observed component 
of the cosmic-ray spectrum beyond the GZK energy. 

The channel in eq. (|2Uj) differs from the previously analyzed resonant selectron and squark 
production under various aspects. First, due to the lightness of the target v mass, it requires 
much larger x\ energies according to eq. JHJ). Second, optimizing cross sections by considering 
scenarios where the sneutrino is maximally coupled to the initial state, and, at the same time, it 
decays into the initial particles v — > \\v is not viable, since this would give rise to invisible final 
states. On the other hand, the sneutrino can decay into heavier neutralinos \\ (* = 2,3,4) 
and charginos xf (j = 1)2). The subsequent neutralino and chargino decays into hadrons 
and leptons might produce highly energetic protons, pions and electrons that would be able 
to initiate the super GZK air showers. Looking back at eq. (jUJ) for the cross section, one is 
interested in scenarios where neither B[y — > Xi u ) nor B(p — > X V i S ) are small. The latter are 
complementary quantities, since one can assume quite accurately B[y — > X V i S ) ~ 1 — B{y — > 
Xi^)- Hence, the sneutrino production cross section is automatically penalized. Furthermore, 
the rates will be in general quite model dependent. Our results will not be simply determined 
by the particle mass spectrum (as in the right selectrons/squarks production), but will depend 
also on couplings (and on the tan/3 and fi parameters) in a nontrivial way. Finally, since a 
visible sneutrino channel will arise from a decay chain, the energy of the visible decay products 
can in principle be quite depleted. In the following, we will also discuss the sneutrino energy 
flow in the visible decay products. 

6 The distances of the center of dense clusters from the earth are estimated to be 16 Mpc for Virgo, while 
our Super-Galactic Plane diameter is about 50 Mpc |40| . 
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Model 


I 


L 


mi/2 (GeV) 


350 


450 


m (GeV) 


180 


303 


tan p 


35 


47 


sign(/i) 


+ 


+ 


OA 2 


n 1 1 7 
U. 11 1 


n no 
U.llo 


mi, (GeV) 


291 


A Ci A 

424 


m x o (GeV ) 


138 


181 


m x o (GeV J 


265 


orn 

350 


m x ± (GeV) 


zoo 


orn 

350 


Y 9 (MeV) 


329 


885 




n 71 

U. 1 1 






n dq 


n i q 


£? x ± (GeV) 


0.20 


0.39 


<jf ak (nb) 


39.4 


19.8 


af a *T 5 (nbGeV) 


13.0 


17.5 



Table 3: Relevant mass spectra, widths, branching fractions and peak cross sections for the 
resonant sneutrino scattering, in the two mSUGRA scenarios I, L defined in Table [0 

The sneutrino branching fraction into visible final states Biy — > X V i S ) is in general made 
up of different components 

B{y - X vis ) l x f) + ££(*>- i/x?) 5(x° - //x?) ~ 1 - fl(i> -> X » (21) 

j=l i=2 

where f ^ v. 

The i> decays into the second neutralino x\ an d into the lightest chargino Xi 1 , when allowed 
by phase space, are in general dominant. The subsequent decays Xi ~^ QQ'Xi (with a branching 
ratio of about 66 % ) and Xi ~~ > l^iXi (with a branching ratio of about 33 % ) are the sources 
of hadrons and charged leptons, respectively, that may trigger the highest-energy air showers. 
The second neutralino can have a substantial decay rate into the invisible channel x\ ~^ u ^Xv 
Otherwise, a x° decay produces quarks and leptons via x° ~~ * QQXii MXi- The X2 ~ * e &R decay 
may also be kinematically allowed. In the latter case, no hadron would arise from the decay 
chain. 

The presence of a heavy invisible particle (xT) among the v decay products makes the visible 
part of the spectrum quite different from the Z-burst case. The results are in general model 
dependent and less promising, and the incoming UHE x° energy should be at higher energies. 
The decay kinematics of a Z-burst event is quite simple. The hadrons resulting from the Z — > qq 
decay carry away all the energy of the incident UHE neutrino (i.e., the Z resonance energy). 
On the other hand, the average energy of the visible sneutrino decay products is lower than 
the sneutrino mass, since hadrons can only appear either from a three-body decay of X21 Xii or 
from a two-body decay chain [xlKxi 1 ) ~~ > Z(W)xi, with Z(W) — > hadrons] following the initial 
two-body v decay. Sneutrino masses, anyway, can be quite heavier than the Z mass, and this 
could compensate partly the depletion of the average energy of the sneutrino decay products, 
that in any case would have a less resonant structure than in the Z-burst case. The detailed 
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energy spectrum of the hadrons and leptons coming from the sneutrino decay cascades also 
depends crucially on the masses of the particles in each step of the decay chains. Their study 
would require a Monte Carlo treatment of the complete kinematics. We will not go into this 
analysis here. 

The observed event number is ruled by the product a~, eak To, that is the peak cross section 
(depending on the sneutrino branching fractions into the initial and final states) times the 
resonance width. 
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0.001 
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L °gl0( E jc,v) 

Figure 4: Cross sections (in nb) for resonant sneutrino production and visible decay (solid 
curves) coming from UHE neutralino scattering off relic neutrinos in the two mSUGRA sce- 
narios I, L defined in Tabled The dotted line represents the resonant Z-boson cross section for 
the scattering of UHE neutrinos off relic neutrinos. The dot-dashed line refers to the continuum 
cross section for ViV T ^ l%c — > ijj. We assumed a relic-neutrino mass m„ = 0.1 eV. 

In Table El we present these quantities, along with the relevant mass spectrum and branch- 
ing fractions, for the scenarios I and L, that are the only benchmarks among the ones presented 
in Table ^ for which sneutrino visible decays are allowed by phase space. In particular, sneu- 
trino decay into the lightest chargino and second lightest neutralino (that are almost degenerate 
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in mass) are open. The sneutrino decay into heavier charginos and X3 4 neutralinos are not al- 
lowed, which is important in order not to deplete too much Biv — > Xi 1/ ) an d, as a consequence, 
the sneutrino cross section. 

In Fig. HI we plot the cross sections for the resonant v production in the channel in eq. (f2T)|) 
versus the incident neutralino energy, for the scenarios I and L. The cross section for the Z- 
production from UHE-z/ scattering off relic neutrinos is also shown versus the incident neutrino 
energy, for comparison. A relic neutrino mass of 0.1 eV is assumed in all cases and the 
resonant incoming UHE neutralino energy is given by E^ ak ~ (m~ — m?.)/(2m I/ ). 

Both the sneutrino peak cross section and the sneutrino width in scenarios I and L are 
quite smaller than (although not very far from, cf. Table EJ) the corresponding Z quantities, 
which are a v ^ ak ~ 185 nb and Y z — 2.5 GeV. In particular, one gets <7~ eafc • T„ ~ 13 — 18 nb 
GeV, to be compared with a v z eak ■ Y z — 463 nb GeV. 

Of course, the latter comparison assumes a comparable flux of incident neutrinos and 
neutralinos. With this hypothesis, the Z-burst is expected to be dominant on the sneutrino- 
burst. Nevertheless, an extra fraction of sneutrino events could add to the Z-burst counting, 
presumably with a bit less structured energy distribution. Further in-depth analysis are needed 
to establish the latter point. 

Note that, while calculating the peak cross sections, in both the Z and v cases the relic 
light neutrinos are assumed to be Dirac fermions. In case of Majorana neutrinos, the peak 
cross sections are enhanced by a factor 2, doubling the number of the events in each case. 

7 Conclusions 

If elementary-particle interactions are governed by supersymmetry, then an important compo- 
nent of the UHECR may be given by the lightest stable supersymmetric particle, namely the 
lightest neutralino. For instance, some of the models trying to explain the UHECR spectrum 
beyond the GZK cut-off via super-heavy particle decays, when including supersymmetry, pre- 
dict a flux of UHE neutralinos arriving to the earth. One of the major challenges of the future 
cosmic-ray experiments is to verify such predictions. 

In this paper, we computed the event rates for right selectrons and right squarks produced 
on resonance, when UHE neutralinos scatter off the electrons and quarks in a detector like 
IceCube. We assumed a well-motivated model of supersymmetry, that is minimal supergravity, 
and analyzed the signal for moderate scalar masses, when the resonant scattering is expected 
to be dominant on the t-channel scattering. 

The right scalar sector turns out to be particularly straightforward to analyze. One has in 
general B{e R — > e X\) — 1> an d B(q R — > q X\) — 1- As a consequence, after constraining the 
mSUGRA parameter space by negative collider searches and neutralino relic-density bounds, 
the relevant phenomenology depends only on the physical masses of the resonant right scalar 
and lightest neutralino. Any other parameter dependence, such as the one arising from the 
neutralino physical compositions and couplings, drops off. We performed a detailed study in 
scenarios inspired to the post-WMAP mSUGRA benchmarks of [33J. 

For comparison, event rates for the expected resonant W signal in UHE anti-neutrino 
scattering off electrons in matter were also presented. 

Two different phenomenological forms for the neutralino fluxes were assumed to calculate 
the event rates: power laws with exponent either /3 = 1-5 or (3 = 2. We set the flux 
normalization on the basis of present bounds on the UHE neutrino flux from the AGASA and 
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Fly's Eye experiments (assuming, quite conservatively, that interaction rates of neutrino and 
neutralino with matter are of comparable strength). 

For similar neutrino and neutralino fluxes, in the leptonic channel the event rates for 
SUSY particle resonances are very promising (up to a few tens of events per year in IceCube). 
They are remarkably greater than the leptonic Glashow W signal for most of the post-WMAP 
benchmarks characterized by moderate scalar masses (m e - s < 500 GeV), that fall either into the 
'bulk' or into the co-annihilation region. The corresponding hadronic signal is less structured 
in the energy spectrum, and penalized by the larger resonance masses (rriq R > 600 GeV). A 
few events per year in Icecube are expected in the most promising benchmarks, in the hadronic 
case. Signal in the 'funnel' and 'focus point' benchmarks is depleted by the large m down to 
less than 0.01 event/year. 

We showed here for the first time the detailed spectrum of the visible energy released in 
the v e e~ reaction in underground detectors, and discussed how to calibrate the SUSY signal 
through the expected W spectrum. 

We also considered the possibility of resonant production of sneutrinos in the interaction 
of UHE neutralinos with light relic neutrinos in the extended Galactic or Local Group hot 
neutrino halo. The visible decay products of the sneutrinos might trigger air-showers beyond 
GZK energies, thus mimicking the so called Z- b urst process. However, even in the most promis- 
ing scenarios, we found that the sneutrino-burst event rate would be more than one order of 
magnitude smaller than the Z-burst rate (assuming similar fluxes for UHE neutrinos and neu- 
tralinos). Furthermore, the UHE neutralino energies tuned for the sneutrino-burst axe usually 
higher than for neutrinos in Z-burst, making it less attractive than the standard Z-burst model. 

In conclusion, we found a well-defined and interesting window of mSUGRA parameters that 
can be tested by the next-generation underground detectors. The corresponding signal may 
exceed in some cases the standard neutrino ones. Its imprint may be disentangled in a realistic 
way. These high-energy astrophysical traces could offer a new opportunity to high-energy 
astrophysics to anticipate SUSY discovery at colliders. 
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